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Abstract. A radio frequency (RF) transition is used to convert a pure F = 2, mp = 2 
87 Rb Bose-Einstein condensate confined in a TOP trap to a mixture oi F = 2,itif = 2 
and F = 2, rap — 1 states. We show that the nature of this coupling process is 
strongly influenced by the presence of the time varying field of the TOP trap, and 
complicated by the presence of multiple Zeeman substates. In particular, the effective 
Rabi frequency associated with the applied RF field is not constant across the spatial 
extent of the cloud leading to a complex geometry for atom-laser output coupling and 
'averaging out' of Rabi oscillations. Further a time-varying detuning can give rise to 
complex spatial structures. 
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1. Introduction 

The development of atom-laser output coupling schemes and techniques for coupling 
single state Bose-Einstein condensates (BEC) into two-component condensates have 
been the subject of many recent experimental studies. The motivation behind these 
studies is the development of practical output couplers and the investigation of 
interactions within condensate mixtures. For these experiments, robust and controllable 
production of different output states is essential. Two-component condensates have been 
produced using sympathetic cooling H, two photon (RF and microwave) transitions 0, 
and RF coupling in optical dipole traps |§. The later allows simultaneous trapping of 
all the states within a hyperfine manifold. Several output coupling schemes have been 
demonstrated experimentally: simple RF output coupling ||, ||, |], optically induced 
Bragg diffraction J7|, || and Josephson tunneling |J. Of the six research groups who 
performed these studies, four used time-averaged orbiting potential (TOP) traps [|Kj 
but none considered effects introduced by the time-varying field. 

This paper extends the study of RF output coupling to include dynamics associated 
with the time-varying field of a TOP trap. RF coupling in a static trap was recently 
characterised by jl2| , and a high degree of control demonstrated. In a TOP trap, the 



population dynamics are strongly influenced by the time- varying nature of the trap. We 
find that the presence of the rotating magnetic field can significantly alter the apparent 
spectral width of the applied RF field, which in turn changes the spatial profile of the 
coupling region. We also find that the large quadrupole gradient of the TOP trap gives 
rise to a large change in the effective Rabi frequency across the cloud, leading to complex 
spatial population distributions. 

Our experimental procedure has been described in ||. We start with a 87 Rb Bose- 
Einstein condensate in the \F = 2, nip = 2) state, confined in a TOP trap with radial 
quadrupole field gradient B' q = 200 G cm -1 and bias field -Btop = 4 G. We then apply 
an RF coupling field with polarization perpendicular to the plane of rotation of the 
TOP field. The effect of this field is to induce transitions between adjacent mp states 
in the F = 2 hyperfine manifold, and hence the initially pure |2, 2} condensate state 
becomes a mixture of the five possible mp states. By changing the frequency, duration 
and magnitude of the coupling field different admixtures of the manifold states can 
be generated. Note that states |2, —1) and |2, —2) are high field seeking and therefore 
actively repelled from the trap, and state |2, 0) being essentially magnetically insensitive 
is untrapped; hence these states will be present in the trapping region for only a few 
milliseconds after the coupling pulse. Once the coupling field has been applied we 
separate the trapped components by removing the trapping potential and 1 ms later 
applying a magnetic field gradient which introduces a Stern-Gerlach separation. The 
different magnetic sub-states (|2, 2) and |2, 1)) become easily resolvable after about 5 ms. 
Although the probing pulse is tuned on resonance with the |2,2) state, the population 
of the |2, 1) state is also accurately determined since the 50 /zs probing pulse (in a 4 G 
bias field) leads to rapid optical pumping into the |2,2) state. This gives us effective 
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resonant probing of both condensate mp states. Typical results of this process are shown 
in Fig. ffl. 





Figure 1. (a) An image of a Bose condensate without any coupling field applied, 
(b) A similar condensate after application of an RF field tuned to the centre of the 
trap with magnitude 72 mG applied for 32 ^s. The individual images are 250 /jm 
along each axis and are taken after a 5 ms time-of-flight through the inhomogeneous 
magnetic field. 



2. Coupling in the TOP trap 

To analyze this coupling process, it is helpful conceptually to first consider the case 
of a static magnetic trap, such as a Ioffe-Pritchard (I-P) configuration. When an RF 
field is applied to couple atoms into other states, the static magnetic field gradient gives 
rise to a static but spatially varying detuning A. To calculate the linewidth of a pulse 
of duration r pulsc we must consider the relative effects of the pulse spectral width and 
power broadening (note that the natural linewidth of an RF transition is insignificant 
in comparison). The pulse spectral width is given by 

r ~ 271 

1 spectral ~ j 
Tpulsc 

and the power broadening (for |2,2) <-> |2, 1)) by [IT 

powor " V2h ' 

where _B r f is the magnetic field amplitude of the RF coupling field. For the |2, 1) <-» |2, 0) 
transition the linewidth is a factor of \j2/3 smaller. 

These two mechanisms combine to give an intrinsic linewidth associated with the 
RF pulse. The detail of how these combine depends on the exact form of each lineshape. 
The lineshape associated with power broadening is Lorentzian. However, the spectral 
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lineshape is determined by the shape of the RF pulse, which for a square pulse gives 
rise to a sine spectral distribution. Numerical simulations suggest that addition in 
quadrature is a good approximation, so that 



■Tpulse ~ y r s p ec t r al ~\~ Fp 0wer . 

We now include the effects of a time varying magnetic field. In a typical TOP 
trap significant angular change in the orientation of the bias field will occur during 
even a short RF pulse (e.g. with a bias field rotating at c^top = 7 kHz, a 10 fis pulse 
corresponds to a 25° change in direction). The effect of this rotation is to introduce time 
dependence to the already spatially varying detuning. This time dependence affects the 
coupling process in two distinct ways. Firstly, as the quadrupole field rotates about 
the centre of the cloud the detuning at any fixed point varies sinusoidally with time. 
Secondly, this sinusoidal variation in the detuning introduces a spatially dependent 
transit time broadening of the transition linewidth. Effectively, the available interaction 
time is reduced from r pulse to r transit by the time varying detuning. This transit time 
broadening, r transit, depends on the distance r from the axis of rotation of the TOP field 
(in our case the z-axis) and r pu i se 

Ttransit = 27rr^ TO P Q , (1) 
J- pulse 

where u' q is the radial gradient of the RF transition frequency (arising from the Zeeman 
shift associated with the quadrupole field gradient B' q ) given by 

4 = 9*Kt- (2) 
By careful choice of experimental parameters, one might choose to avoid the 
complications introduced by the rotating field. This could be achieved, for example, 
by broadening the RF transition so that the coupling rate is essentially the same for 
all atoms. In our experiment, the transition frequency gradient is 14 kHz /im -1 and 
our condensate diameter is ~ 30 /im, giving a change in the transition frequency across 
the spatial extent of the cloud of 420 kHz. Therefore, in order for the coupling to 
be insensitive to the field rotation, the coupling field would need to have a linewidth 
significantly larger than 420 kHz. This could be achieved in two ways, with a very short 
pulse (duration <C 2/is) or a very intense RF field (amplitude 3> 0.85 G). 
Using Eq. (H) we can now determine the total RF transition linewidth 



r ~ y Fp U i se -|- rt rans it. 

The spatial dependence of the linewidth is shown in Fig. for three different values of 

Tpulse- 

There are two regimes that simplify the expression for T. The first of these is where 
Ttransit <S Tpuise, m which case transit broadening can be ignored, and the linewidth of 
the coupling process can be treated as a constant (e.g. Figure 0, r pu i se = 8 /i s). 
The second case is when r transit r pulse (e.g. in the case of quasi-continuous output 
coupling). From Eq. ([!]) we see that the linewidth in this case is approximately linear 
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Figure 2. The total linewidth T as a function of position within the condensate for 
weak RF field and three different pulse lengths: — 32 /j,s, 16 (is, and • • • 8 /is. 



with radius. However, in our experiments neither of these simple approximations apply 
and a more complete treatment is required. 



3. Modeling the coupling 

We are able to broadly predict the population dynamics of the various mp states as a 
function of time with the following simple Rabi cycling model that neglects any kinetic 
or mean-field effects within the condensate. In an F = 2 manifold there are four Rabi 
frequencies Q1...4 associated with the RF field, where £l\ couples states |2, 2) <-> |2,1), 
Q2 states |2, 1) <-> |2,0), etc. From the Clebsch-Gordan coefficients for the transitions, 
Qi = fi 4 , n 2 = 1^3 and v^3^i = \f2Vt 2 allowing us to define a single Vt{t) for all 
transitions 

£F/i B -B r f(t) 



tt{t) 



V2h 



Each RF transition will be detuned an amount A(r, 6, t) = 8(r) cos(o;top£ + #) from 
resonance. The time dependence in A is introduced to account for the rotation of the 
TOP field, and 5(r) accounts for the effects of uj' q introduced in Eq. (@). Using these 
parameters we can write the following equation for the population amplitudes of the 
m F states a 2 ..._2 



d_ 

dt 



< a 2 \ 

a 1 
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a-2 / 
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tt(t)* 2A(r,9,t) 
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This model ignores trap loss, this is because we estimate that the average time 
taken for an untrapped atom to be ejected from the trap is approximately 3 ms, so the 
loss rate 7 (~ 0.3 kHz) is small compared with typical values for Q and A (in the range 
10 to 100 kHz). 
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Time ((Is) Time ((Is) Time ((Is) 

Figure 3. The populations of the |2, 2) ( ) and |2, 1) ( ) states as a function 

of time for Q = 50 kHz : (a) homogeneous coupling, (b) inhomogeneous coupling with 
very slow bias field rotation and (c) inhomogeneous coupling with bias field rotation. 

Figure |3] shows the calculated populations of states |2, 2) and |2, 1) as a function 
of time for three distinct regimes: (a) strong coupling such that the coupling rate is 
uniform across the condensate, (b) very slow bias field rotation and weak (and therefore 
inhomogeneous) coupling across the condensate, and (c) typical experimental conditions 
of weak coupling with rapid bias field rotation, which introduces temporal variation 
into the coupling. Unless stated otherwise, to match our experimental conditions we set 
S(r) = o; q r with uj' q = 14 kHz /mi -1 , and allow a 1 fxs rise and fall time on B T f(t) which has 
an amplitude of 72 mG (equivalent to Q ~ 50 kHz). To calculate the total population in 
each state we assume a 3-D Thomas- Fermi distribution of size 30 fim x 30 /im x 11 /j,m 
and integrate over this region. 

When the linewidth is significantly larger than the change in transition frequency 
across the condensate (as in Fig. ||(a)), spatial effects are reduced and the population 
essentially Rabi cycles between the fully stretched states |2, 2) and |2, — 2). For ease 
of comparison we keep Q constant across all plots in Fig. [3] and simulate the broad 
linewidth condition in (a) by reducing the quadrupole field gradient by a factor of 200. 
For Figs. ^|(b) and (c) we return to conditions more relevant to our experiment. For Fig. 
|3](b), although we have eliminated bias field rotation, complex spatial behaviour results 
since the coupling linewidth is smaller than the Zeeman shift across the condensate. 
However the decay of Rabi cycling due to large detuning changes across the cloud leads 
to an equilibrium population distribution that depends only on the pulse linewidth of the 
transition relative to the Zeeman shift across the condensate. Figure |^(c) has identical 
parameters to Fig. |3](b) except that the rotation of the bias field has been included, so 
that the time dependence of A is significant. In this case the Rabi cycling is rapidly 
overwhelmed by a population oscillation at the bias field rotation frequency. This cycling 
appears only when r u' q ^> ^> cutop where r is the size of the condensate, and can 
be thought of as a frequency modulation of the coupling field with a beat note at the 
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Figure 4. Spatial structure of the |2, 2) state after RF coupling for (a) 6 {is and 
(b) 136 /is. The condensates shown are 30 /an across. 

bias field rotation frequency, leading to slow cycling of the population. Along with this 
cyclic behaviour, our modeling shows that the coupling also introduces interesting small 
scale spatial structure in the condensate such as shown in Fig. f|. 

4. Experimental results 

The focus of this experimental work is to understand the spatio-temporal dynamics of 
RF coupling in a TOP trap. To test the behaviour predicted by the model we first 
produce Bose condensates in the |2, 2) state as described earlier, and then by varying 
the duration and power of the RF coupling field we hoped to obtain results similar to 
those in Fig. |^(b) and (c). 

We first considered the case of a relatively short RF pulse, since this minimises 
the effect of bias field rotation on the coupling surface. To do this we investigated the 
relative populations of the |2, 2) and |2, 1) states as a function of RF power for a fixed 
pulse length of 20 /is (compared with a trap rotation period of 143 /is). The results of 
this, along with a simulation, are shown in Fig. |5|. It should be noted that there are no 
free parameters in any of the simulations of relative population. 

The agreement between the experimental and theoretical results is surprisingly good 
given the simplistic nature of the model. At a field amplitude of 150 mG we obtain an 
approximately equal proportion of atoms in the two trapped states. For comparison we 
note that this field is significantly larger than that typically required for evaporative 
cooling. Within the limits of experimental uncertainty, the decrease of the |2, 2) state 
population with increased RF field behaves as we might expect, whilst the |2, 1) state 
population appears to grow more rapidly than predicted. The reason for this minor 
discrepancy may be due to the quadratic Zeeman effect which would inhibit cycling into 
the untrapped states. 
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Figure 5. Populations of the \2, 2) (+) and |2, 1) (o) states as a function of the applied 
RF magnetic field. The lines are the results of simulations described in the text. The 
error in the data is of order ±10% of N . 



Secondly, we considered the case of a weak pulse, since this reduces the effect of 
power broadening on the linewidth. To do this we investigated the relative populations 
of the two trapped states as a function of RF pulse length for an RF field amplitude 
of 72 mG. We varied the pulse length from 10 /zs to 1 ms to cover a range of times 
relative to the bias field rotation period. Results are shown in Fig. |6|. The parameters 
match those used to generate Fig. §(c). Rabi cycling is suppressed as expected, and 
the relative populations rapidly become constant and equal. However, the data does 
not exhibit the cyclic behaviour of the simulation and we also observe a linear increase 
in internal energy of the condensate as a function of RF pulse length. We believe that 
both these observations stem from the spatial structure introduced by the coupling 
field. There are large phase gradients associated with the predicted spatial structure 
and these correspond to increased internal energy. The simulation does not include 
diffusion effects, and atomic movement on the scale of the spatial structure will tend to 
average out the slow population cycling at the bias field rotation frequency. 

5. Conclusions 

To conclude, we have characterised the RF coupling process for Bose condensates in 
a TOP trap and shown that it is possible to engineer predetermined multi-state Bose 
condensate mixtures. A simple model predicts the essential features of the condensate 
behaviour. Our simulations indicate that the coupling process introduces complex 
spatial features. Preliminary results using a full GPE calculation indicate that simple 
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Figure 6. Populations of the |2,2) (+) and |2,1) (o) states as a function of pulse 
length. The error in the data is of order ±10%. 



RF coupling in a TOP trap can lead to vortex formation fllci |. This will be the subject 
of future work at Otago. 
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